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C hern- Simons Vortices 



C. DUVAL (1) and P. A. HORVATHY (2) 



Abstract. In (2 + 1) dimensions, the Maxwell term —jFapF'^^ can be replaced by the Chern-Simons three- 
form {K/A)e"^''AaFi3^, yielding a novel type of ' electromagnetism' . This has been proposed for studying 
the Quantum Hall Effect as well as High- Temperature Superconductivity. The gauge field can be coupled 
to a scalar field either relativistically or non-relativistically. In both cases, one admits finite- energy, vortex 
' solutions. 

o ■ 

^ ■ Contemporary Mathematics 203, pp. 271-288 (1997) 



, 1. Introduction 

The phenomenological description of 'ordinary' superconductivity is provided by Landau-Ginzburg the- 
^ ' ory [1]: the Cooper pairs formed by the electrons are represented by a scalar field, whose charge is twice 

■ that of the electron. The scalar fields interact through their electromagnetic fields, governed by the Maxwell 

. equations. The theory admits static, vortex-like solutions [2]. 

, Landau-Ginzburg theory is non-relativistic; a relativistic version is provided by the Abelian Higgs model 

fT^ I [3], which admits again static and purely magnetic vortex- type solutions. 

In the eighties, however, new phenomena were discovered, namely the Quantum Hall Effect [4] and 
High- Temperature Superconductivity [5]. None of these is described by Landau-Ginzburg theory. Hence the 
need for a new theory. 



o 



Oh! 

, Both phenomena mentioned above are essentially planar. The mathematical room for a modification is 



provided by the three-form 



Xl : (1^1) je"^A.Fg„ 

introduced by Chern and Simons [6] in their study of secondary characteristic classes. This form, which only 
exists in (2 + 1) dimensions, can be added to the conventional Maxwell term —jFa/sF"^ [7]. The use of the 
Chern - Simons term has in fact been proposed in the above contexts [8] , [9] . 

Since the Chern - Simons term is of lower-order in derivatives, for large distances it dominates the 
Maxwell term. Thus, for describing long-range phenomena, the Maxwell term can even be dropped: the 
dynamics of the 'electromagnetic' field will be described by the Chern-Simons term (1.1) alone. 

Below, we review various aspects of Chern-Simons gauge theory. First, we explain the construction of 
relativistic topological vortices. The non-relativistic limit is physically relevant owing to the intrinsically 
non-relativistic character of condensed matter physics. It also provides an explicitly solvable model. Finally, 
we consider spinorial models. Again, explicit solutions are found. 



2. Relativistic Chern-Simons vortices 
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Let us consider a complex scalar field, ij), coupled to an abclian gauge field, described by the vector poten- 
tial [1-form], Afj^dx^, defined on (2-|-l)-dimensional Minkowski space with the metric (g^i/) — diag (1, — 1, — 1), 
the coordinates being =t and {xi) = x. The action S = J d^xC is given by [10] 

(2.1) C = \D^^{D^i^Y + \n e"^^A„J>^ - C/(^), 

where -D^n'^i' = 9(nV ~ i^A^ip denotes the covariant derivative, the constant e is the electric charge of the field 
ip, and K is a constant. The potential f/(V') is chosen, for later convenience, to be 

(2.2) u{^P) = ^\^Pm\^P\'-l)\ 

The Euler-Lagrange equations of the action (2.1) read 

SU A 

^V^U-ip = - 

(2.3) 



where Fa^ = daAjj — d^Aa is the 'electromagnetic' field, and = {q,j) is the current 

(2.4) j/^=|[^*£,''V-^(^^V^)*]. 

The first of the equations in (2.3) is the Non-Linear KIcin-Gordon equation (NLKG) familiar from the 
Abelian Higgs model [3]; the second, called Field-Curent Identity (FCI), replaces the Maxwell equations. 

It follows from the Bianchi identity, e"^"*"' daFpj = 0, that the current (2.4) is conserved, <9^j^ = 0. 

Due to the presence of the vector potential in the Chern-Simons term, the Lagrangian (2.1) is manifestly 
not invariant with respect to a local gauge transformation 

(2.5) ^{x) ^ e^^^("V(^), A^^A^ + d^X. 

Using the Bianchi identity, the change can be written, however, as da (('«/4) e"^'' Fp^ A) . This is a surface 
term, which does not change the equations of motion. The two other terms in (2.1) are gauge-invariant. 

The same conclusion can also be reached by looking at the field equations: the FCIs are invariant, since 
both the field strength and the current are invariant; both sides of the NLKG get in turn be multiplied by 



Finite-energy configurations 

Let us consider a static field configuration (A^jV')- The energy, defined as the space integral of the 
time-time component of the energy-momentum tensor associated with the Lagrangian (2.1), is 



(2.6) 



E \ 



J d^xT'^^ = J d^x[^DiiP{D'iP) 



2^ ^ol 



+ kAoB + U{i}) 



where B = —F^^ is the magnetic field. Note that this expression is not positive definite. Observe, however, 
that the static solutions of the equations of movement (2.3) are stationary points of the energy. Variation 
of (2.6) w. r. t. Aq yields one of the equations of motion, namely 



(2.7) 
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Eliminating Aq from (2.6) using this equation, wc obtain the positive definite energy functional 
(2.8) E = jd?r 



We are interested in static, finite-energy configurations. Finite energy at infinity is guaranteed by the 
conditions {}) 

= o(l/r), 
= o(l/r), r ^ oo. 











(2.9) 


< 


a.) 


B 






^ Hi.) 





= o(l/r). 

Therefore, the U{1) gauge symmetry is broken for large r. In particular, the scalar field is covariantly 
constant, Dtf) = 0. This equation is solved by parallel transport, 



V'(a;) = exp 



I I eAidx^ 

So 



^0, 



which is well-defined whenever 



(2.10) 



eAidx'' 



<fxeB = e^ = 2'Kn, n = 0, ±1,... 



Thus, the magnetic flux is quantized. 

By i.), the asymptotic values of the Higgs field provide us with a mapping from the circle at infinity into 
the vacuum manifold, which is again a circle, = 1. Since the vector potential behaves asymptotically as 



- --djlogtp, 



(2.11) 

the integer n in Eq. (2.10) is precisely the winding number oi this mapping; n is called the topological charge 

(or vortex number). 

Spontaneous symmetry breaking generates mass [12]. This occurs in a somewhat unusual guise, though. 
Expanding around the vacuum expectation value of tp we find = —eA^ so that the FCI in (2.3) is 
approximately 

Hence Fap. = — (e^/k) e^f^pA^ . Inserting here Fa/s and deriving by d" we find that the gauge field A^ satisfies 
the Klein- Gordon equation 

(2.12) DA" = j A", 
showing that the mass of the gauge field is 

(2.13) ruA = — . 

The Higgs mass is found in turn in the usual way: the expectation value of the Higgs field can be chosen 
as Vo = (Ij 0). Expanding ip as {tpr, V*??) = (! + '/') ^) we get 



(2.14) 



since 



(2.15) 



=0 ^ 



1 ffU_ 



=0 



1 is a critical point of U. We conclude that the mass of the Higgs particle is 



IVI=1 



2A. 



(This result can also be seen by considering the radial equation (2.17) below). 



(^) These conditions are in no way necessary; they yield the so-called topological solitons. Non-topological 
solutions are constructed in Ref. [11]. 
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Raclially syiiiiiiotric solutions [10] 

For the radially symmetric Ansatz 

(2.16) A^^A^{r), A, = Q, = A{r), ^{r) = J {r)e-'^\ 

the equations of motion (2.3) read 

r ar k 
with asymptotic conditions 

Tl 

liuir^oo A{r) = , limr^oo/(r-) = 1, 

(2.18) e 

limr^o A{r) = 0, hmr^o/('') = 0. 



This is either seen by a direct substitution into the equations, or by re- writing the energy as 



1r 2e4/2 

where a = eA + n. The upper equation in (2.17) is plainly the radial form of (2.7). Then variation of (2.19) 
with respect to a and / yields the two other equations in (2.17). 

Approximate solutions can be obtained by inserting the asymptotic value, / ~ 1, into the first two 
equations: 

«' , ^ </ a „ 

- + —Ao = 0, A'o + =0, 

r K K r 

from which we infer that 

(fAo 1 dAo _ / 2 / ^ 

dp-' p dp 

This is the modified Bessel equation [Bessel equation of imaginary argument] of order zero. Hence 

(2.20) Ao = CKoi-r). 
Similarly, fov a = n/e + A we find, putting a = a/r, 

a" + — - fl + ^)a = 0, 
p V p-'/ 

which is Bessel's equation of order 1 with imaginary argument. Thus a = CK\{p) so that 

(2.21) A = --+C-rKA-r). 

Another way of deriving this result is to express A from from the middle equation in (2.17), 



n K ^ ^ 
e dr ^ 
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The consistency with (2.21) follows from the recursion relation Kq = —Ki of the Bcsscl functions. 

An even coarser approximation is obtained by eliminating the ^ term by setting a = ut~^I'^ and 
dropping the terms with inverse powers of r. Then both equations reduce to 



(2.22) 



U = — U 
^ K ' 



c _ 

\q = a = —=e 



which shows that the fields approach their asymptotic values exponentially, with characteristic length deter- 
mined by the gauge field mass. 

The deviation of / from its asymptotic value, y> = 1 — /, is found by inserting into the last eqn. of 
(2.17); developping to first order in if> we get 



(2.23) 



1 



^p<' + _(^' _ 2A<^ ~ 



(p = C7i<:o(V2Ar), 



whose asymptotic behaviour is again exponential with characteristic length (to,/,) ^, 



(2.24) 



The penetration depths of the gauge and scalar fields are therefore 



(2.25) 



1 ^ ^ , ^ 
7] = = — and t, = 

niA K 



1 



2A 



respectively. For small r instead, inserting the developments in powers of r, we find 

fir) ~ /orl"l + . . . , 



(2.26) 



^o~ao-#^r2|"l+... 

2k n 



e /qQ^O ^2|n|+2 

2K(|n| + 1) 



+ 



where ao and /o are constants. In summary, 

IV'(OI = /W oc 

(2.27) \E{r)\ = \A',{r)\ <x 

\B{r) 



r 



oc 



'rl"l r ^ 

fr^l"!-! r 

[ Cr~^/^ e"'"^'' + lower order terms r — > OO 

r^l"! r ~ 

^-ruAT _|_ j^^gj. r ^ 00 



Self-dualitv 

Let us suppose that the fields have equal masses, = m\ = and hence equal penetration depths. 
Then the Bogomolny trick applies, i.e., the energy can be rewritten in the form 



(2.28) 



I 



\\{Di±iD2) 
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The last term can also be presented as 

The integrand of the _B-tcrm yields the magnetic flux; the second is transformed, by Stokes' theorem, into 
the circulation of the current at infinity which vanishes since all fields drop off at infinity by assumption. Its 
integral is therefore proportional to the magnetic flux, ±e$/2. Since the flrst integral is non-negative, we 
have, in conclusion, 

(2.29) E>^=7r\n\, 
equality being only attained if the 'self- duality' equations 

2 

(2.30) Ditp = TiD2i^ and = ±^ |^|2(1 - |^|2) 

hold. It is readily verified that the solutions of equations (2.7) and (2.31) solve automatically the non-linear 
Klein-Gordon equation. 

For the radial Ansatz (2.16) the self-duality equations become 



(2.31) f' = ±-f, - = ±X/'(/'-l)' 

r r 



where we introduced again a = eA + n. Deriving the first of these equations and using the second one, we 
for / we get the 'Liouville - type' equation 

2 

(2.32) Alog/ = ^/2(/2-l). 

Another way of obtaining the first-order eqns. (2.32) is to rewrite, for 

[/(/) = (mV8)/2(/2-l)^ 



the energy as 







a „t2 1 rfl' rn^ 



oo 

=F ira 



The boundary conditions read 

(2 33) ""^"^^ 

a(0) =n, /(O) =0, 

and thus E > Tr\n\ as before, with equality attained iff the equations (2.32) hold. 
For n = the only solution is the vacuum, 

(2.34) / = 1, A=l. 

To see this, note that the boundary conditions at infinity are /(oo) = 1 and ^(oo) = 0. Let now /(r), A{r) 
denote an arbitrary finite-energy configuration and consider 



(2.35) 



Mr)=f{r), Ar{r)=TA{r) 
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where r > is a real parameter. This provides us with a l-parameter family configurations with finite energy 



(2.36) 



2r-' 2m^^rf' 



U{f) , 



which is a monotonic function of r, whose minimum is at r = i.e. for a = 0. Then Eq. (2.32) implies that 
/' = so that / = 1 is the only possibility. 

Let us assume henceforth that n^O. No analytic solution has been found so far. To study the large-r 
behaviour, put ip= 1 — f. Inserting / 1, Eq. (2.32) reduces to 



(2.37) 

Prom here we get 



Thus, for large r. 



, a a 2 
ip = — = Tm ip. 

r r 



(p" + -if' - m^if = 
r 

a" a' — m^a = 

r 



<fi = CKo{mr). 
a = CmrKi {mr) 



(2.38) 



/ ~ 1 — C-ftro(TOr) and A~ 1 CmrKiimr) 

e 



with some constant C. For small r instead, Eq. (2.32), yields, to ©(r^l"'"*"^). 



(2.39) 



2(2n + 2)2 



^3|n|+2 _^ 0(r^l"l+^). 



A = 



2(2|n| +2)e 



J2\n\+2 , 



„in+2 



2(4|n|+2)e 



+ 0(r*l"l+^). 



The result is consistent with (2.26) since the constant ao is now ao = m/2e = e/2K. 




Fig. 1. The scalar field of the radially symmetric charge-2 relativistic vortex. 
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Fig. 2. The magnetic field of the radially symmetric charge-2 relativistic vortex. Note that B = where the 
scalar field vanishes: the vortex has a doughnut-like shape. 

Let us mention that Eq. (2.33) is actually valid in full generality, without the assumption of radial 
symmetry. Expressing in fact the vector-potential from {Di zb iD2)ip = as 

(2.40) el = V(Arg^) zb V x log |^| 

and inserting into the second equation in (2.31), we get again (2.33), with l-^l replacing /. Index-theoretical 
calculations show that, for topological charge n, Eqn. (2.33) admits a 2\n\ parameter family of solutions 
[10]. 




3. Non-relativistic vortices 

The non-relativistic limit of the above system is found [13,14] by setting 

where ^ and ^ denote the particles and antiparticles, respectively. Inserting this into the action, dropping 
the oscillating terms and only keeping those of order 1/c, we see that both the particles and antiparticles 
are conserved. We can therefore consistently set ^ = 0. The remaining matter Lagrangian reads 

(3.1) Chatter =t^*D,^i>-^-^ + ^i^^^'f, 

where A = e^/mc|K|. One can show that the theory is indeed non-relativistic [15], [16]. In what follows, we 
put c = 1 . 

Let the constant A be arbitrary. Variation of / Cmatter w. r. t. yields the gauged non-linear 
Schrddinger equation 



(3.2) idt^ 



£)2 

eAo - A**^' 

2m 



A fully self-consistent system is obtained by adding the matter action to the Chern-Simons action (1.1); the 
variational equations of this latter are clearly the FCIs in Eq. (2.3), written in non-relativistic notations as 

(3.3) B = e'^diA3 = -- g, E' = -d,A° - dtA' = -t'^ , 

K K 
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where the current is now 

(3.4) = ie, J) = ^ [^*D^ - ). 

This is legitimate, since the Chern-Simons term is also non-relativistic: it is in fact invariant with respect 
to any coordinate transformation. 

Self-dual vortex solutions 

We would again like to find static soliton solutions. The construction of an energy-momentum tensor 
is now more subtle, because the theory is non-relativistic. This is nevertheless possible [15], [16]. One gets 
the energy functional 

(3.) 

Applying once more the Bogomolny trick, this is also written as 
Thus, for the specific value 

(3.6) A = 



e2 



m\K\ 

of the non-linearity, — the same as the one we obtained above by taking the non-relativistic limit of the 
self-dual relativistic theory — the second term vanishes. Then the energy is positive definite, whose absolute 
minimum — zero — is attained for "self-dual" fields i.e. for such that 

(3.7) {Di ± i£»2)* = 0. 

Separating the phase as \1/ = e^'^'^y^, we get 

(3.8) 1= Vw + X log£., 
where g solves the Liouville equation, 

2e2 

(3.9) Alog^ = ± g, 

K 

cf. (2.33). Physically admissible solutions arise when the r.h.s. is negative. Hence the upper sign has to be 
chosen when ac < and the lower sign when k > 0. Then the general solution reads 

(3.10) .-^I'^l l^'l^ 



e2 [1 + l/IT' 
where / is a meromorphic function. In particular, for 

n 

(3.11) /w = E 



i=l ' - 



where the q and the Zi arc 2n complex numbers, we get a 4n-parameter family of physically admissible 
solutions, with vortex number n. 
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Note that g vanishes at the poles of /; these points represent the 'positions' of the vortices. The 
singularity in g is precisely compensated by the singularity in the phase, leaving us with a regular magnetic 
field [15]. 

The number of independent parameters is in fact 4n — 1, due to a global gauge freedom; index-theory 
calculations [17] show that this is indeed the maximal number of independent solutions. 

For / = c/z" in particular, we obtain the radially symmetric solution 



(3.12) 



g{r) 



4n^|K| 



whose flux is evenly quantized. 



(3.13) 



$ = — (sgK) — X n. 
e 



A few non-relativistic vortices are shown on FIGs. 3-6 below. 




Fig. 3 The non-relativistic radially symmetric N = 1 vortex has a maximum at r = 0. 




Fig. 4 The non-relativistic radially symmetric N = 2 vortex has a 'doughnut-like' shape. For N > 2, the 
particle density vanishes at r = 0. 
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Fig. 5 The non-relativistic N — 2 vortex representing two separated 1-vortices with 'positions' at x = ±1 
and y = 0. 




Fig. 6 The non-relativistic N = A vortex representing two separated 2-vortices with 'positions' at x = ±1 
and y = 0. 

Let us mention that the only time-dependent solutions found so far are either obtained from the static 
solution by the action of the Schrodinger group (the fundamental symmetry group of the theory [13, 14, 15, 
16]), or are stationary [18]. 

4. Spinor vortices 



A. Relativistic spinor vortices 

In Ref. [19] Cho et al. obtain, by dimensional reduction from Minkowski space, a (2 + l)-dimcnsional 
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system. After some notational changes, their equations read 
(4.1) 

{ic^±Da -m)ip±= 0, 
where the two sets of Dirac matrices are 

(4.2) (7g) = (±(l/c)a^^a^-^al), 

and the V± denote the chiral components, defined as eigenvectors of the chirahty operator 

(4.3) ^=(~T 

Observe that, although the Dirac equations are decoupled, the chiral components are nevertheless coupled 
through the Chern-Simons equation. Stationary solutions, representing purely magnetic vortices, are readily 
found [19]. It is particularly interesting to construct static solutions. For Aq = and dfAi = 0, setting 



(4.4) = e' 



-imt 



the relativistic system (4.1) becomes, for c = 1, 

Ke'^diAj = -e{\F+\^ + \G-\''), 

(4.5) 

{Di + iD2) F±=0, {Di- iD2) G± = 0. 

Now, for F± = or G± = 0, these equations are identical to those which describe the non-relativistic, 
self-dual vortices of Jackiw and Pi [15] ! 

B. Non-relativistic spinor vortices 

Non-relativistic spinor vortices can also be constructed along the same lines [20] . Following Lcvy-Lcblond 

[21], a non-relativistic spin ^ field ^ ~ ( ) ^^^^^ ^ ^^'^ X s-re two-component 'Pauli' spinors, is described 
by the 2 -|- 1 dimensional equations 

( (a-D)^ + 2m X =0, 



(4.6) 



Dt^ + i{a-D)x = 0. 



These spinors are coupled to the Chern-Simons gauge field through the mass (or particle) density, 
g= ]$]^, as well as through the spatial components of the current, 

(4.7) J = i{^^ax-X^^^), 

according to the Chern-Simons FCI (2.3). The chirality operator is still given by Eq. (4.3) and is still 
conserved. Observe that $ and x in Eq. (4.6) are not the chiral components of -0; these latter are defined 

by i(i±ir)^± = ±V±. 

It is easy to see that Eq. (4.6) splits into two uncoupled systems for ^+ and Each of the chiral 

components separately describe (in general different) physical phenomena in 2 + 1 dimensions. For the ease 
of presentation, we keep, nevertheless, all four components of ip. 
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Now the current can be written in the form: 
(4.8) J = ^ ($t5^ _ (5$)t$) + V X $t,3$) . 

Using {D ■ BY = + eBa^, we find that the component-spinors satisfy 

= -^\D^ + eBa3\x-^{^-E)^. 



(4.9) 



Thus, $ solves a 'Pauli equation', while x couples through the Darwin term, a ■ E. Expressing E and B 
through the FCI, (2.3) and inserting into our equations, we get finally 



(4.10) 



iDtX = 



2m 2mK 



2m 2mK 



{ax 



If the chirality of ip is restricted to +1 (or —1), this system describes non-relativistic spin +i (~^) fields 
interacting with a Chern-Simons gauge field. Leaving the chirality of ip unspecified, it describes two spinor 

fields of spin ± i, interacting with each other and the Chern-Simons gauge field. 

Since the lower component is simply x = — (l/2m)(i7 • -D)$, it is enough to solve the ^-equation. For 



(4.11) 



and = 



respectively — which amounts to working with the ± chirality components — the 'Pauli' equation in (4.10) 
reduces to 



(4.12) 



Pi 
2m 



± A(*]t*±) 



A = 



2mK 



which again (3.2), but with non-linearity ±A, halfoi the special value A in (3.6) used by Jackiw and Pi. For 
this reason, our solutions (presented below) will be purely magnetic, {At = 0), unlike in the case studied by 
Jackiw and Pi. In detail, let us consider the static system 



(4.13) 

and try the first-order Ansatz 
(4.14) 



— {D^ + eBa3)-eAt $ = 0, 
2m 



J 



--V X At, 
e 



, kB = —eg. 



(£»i ±i£>2)$ = 0. 



Eq. (4.14) makes it possible to replace = Df + D2 by T^B, then the first equation in (4.13) can be 
written as 



(4.15) 



-^eB(Tl + (73) - eAt 
2m 



$ = 0, 
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while the current is 

Now, due to the presence of fj^, both Eq. (4.16) and the second equation in (4.13) can be solved with 
a zero At and J: by choosing <i> = <3>_|_ (<& = $_) for the upper (lower) cases respectively makes (=Fl + (73)$ 
vanish. (It is readily seen from Eq. (4.15) that any solution has a definite chirality). 

The remaining task is to solve the first-order conditions 

(4.17) {Di + iD2)^+ = 0, or (£>i - i£>2)*- = 0, 
which is done in the same way as in Section 3: 

- 1 ^ 2e^ 

(4.18) A = ±—VxlogQ, Alog£> = ± Q. 

A normalizable solution is obtained for when «; < 0, and for ^_ when k> 0. (These correspond precisely 
to having an attractive non-linearity in Eq. (4.12)). The lower components vanish in both cases, as seen 
from 

(4.19) ^ = __L(5^.5)$. 

Both solutions only involve one of the 2 + 1 dimensional spinor fields ip±, depending on the sign of k. 

The physical properties such as symmetries and conserved quantities can be studied by noting that our 
equations are in fact obtained by variation of the 2 + 1-dimensional action given in [20] , which can also be 
used to show that the coupled Levy-Leblond — Chern-Simons system is, just like its scalar counterpart, 
Schrodinger symmetric [15]. 

A conserved energy- momentum tensor can be constructed and used to derive conserved quantities [20] . 
One finds that the 'particle number' N determines the actual values of all the conserved charges: for the 
radially symmetric solution, e.g., the magnetic flux, —cN/k, and the mass, A4 = mN, are the same as for 
the scalar soliton of [15]. The total angular momentum, however, can be shown to be / = half of 

the corresponding value for the scalar soliton. As a consequence of self-duality, our solutions have vanishing 
energy, just like the ones of Ref. [15]. 



(4.16) 



2m 



-V X 



C. The non-relativistic limit 
Setting 

(4.20) ^+ = e-^'"^'*(^J+J and V- = e 

Eq. (4.1) become 

(4.21) 



-ca-Dx = 0, 
iDtX + CCT • Z)i> + 2mc^x = 0, 



where $ = and x~\~ .In the non-relativistic limit mc x » iDtX, so that this latter can be 



dropped from the second equation. Redefining x as x = cx yields precisely our Eq. (4.6). This also explains, 
why one gets the same (namely the Liouville) equation both in the relativistic and the non-relativistic cases: 
for static and purely magnetic fields, the terms containing Di are automatically zero. 

In this paper, we only reviewed the Abelian Chern-Simons theories. Non-Abelian generalizations are 
studied in Refs. [14], [22] and [23]. 
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